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Abstract 

Theoretical model of the coherent control of nuclear spin isomers by microwave radiation has 
been developed. Model accounts the M-degeneracy of molecular states and molecular center-of- 
mass motion. The model has been applied to the 13 CH3F molecules. Microwave radiation excites 
the para state (J=11,K=1) which is mixed by the nuclear spin-spin interaction with the ortho state 
(9,3). Dependencies of the isomer enrichment and conversion rates on the radiation frequency have 
been calculated. Both spectra consist of two resonances situated at the centers of allowed and 
forbidden (by nuclear spin) transitions in the molecule. Larger enrichment, up to 7%, can be 
produced by strong radiation resonant to the forbidden transition. The spin conversion rate can 
be increased by 2 orders of magnitude at this resonance. 

PACS numbers: 32.80.Bx, 33.50.2j, 03.65. 2w 
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I. INTRODUCTION 



Nuclear spin isomers of symmetrical molecules are fascinating objects 0. Their properties 
are determined by nuclei' spins deeply hidden in the molecule. Most known are the hydrogen 
isomers that demonstrate anomalous stability, 1 year at ambient temperature and pressure 
0. Latest results on hydrogen isomers can be found in |3|, ^] and references therein. Many 
other molecules have spin isomers too. But so far their physical properties remain almost 
unknown. This makes investigations of spin isomers fundamentally important. Spin isomers 
have also practical applications, e.g., as spin labels, in isomer selective chemical reactions 
H H , or in isomer enhanced NMR technique |7], || . These applications are developed solely 
with hydrogen isomers. Extension to other molecules needs efficient methods of isomer 
enrichment. For a long time enrichment of only hydrogen isomers was possible. Recently a 
few separation methods for polyatomic molecules have been developed (see the review |5|) 
that has advanced the field significantly. Further progress needs new enrichment methods. 

New approach to the problem of isomer enrichment is based on the use of strong elec- 



tromagnetic radiation. There are two modifications of the method. In the first one [TCJ, [11 
(earlier discussion of the radiation effects see in fll2l), radiation populates specific states of 
one spin isomer situated in the vicinity of states of other isomer that consequently results 



in the enrichment. In the second method [13|, enrichment appears due to combined action 



of population transfer, dynamical shift of molecular levels and light induced coherence. The 
latter method (coherent control of spin isomers) promises to be more efficient. To avoid 
any confusion, we note that the light-induced enrichment resulting from stimulating con- 
version of spin species differs radically from the previously known separation methods, e.g., 
light-induced drift method which separates physically the isomers |9|]. 

There are no proofs yet that light-induced enrichment of spin isomers is feasible. We 
are aware of only one attempt to verify the proposals. It was done by microwave excita- 
tion of rotational transition in CH3F. The experiment did not give a positive result. At 
the time when this experiment was performed, only a qualitative theoretical model of the 
light-induced enrichment was available |TIJ. Presently, the understanding of the underly- 
ing physics has been improved substantially. In view of further experiments in this area it 
is desirable to examine the microwave induced enrichment of CH 3 F spin isomers in more 
detail. This is the goal of the present paper. Existing theoretical models of coherent con- 
trol cannot be used for quantitative analysis directly because of their lacking to account the 
M-degeneracy of molecular states. Account of such degeneracy is another goal of this paper. 



II. QUALITATIVE PICTURE AND KINETIC EQUATION 

Previous analysis has shown that significant enrichment can be obtained if radiation 
interacts with the states having large difference in populations. In this context, microwave 
excitation is not the best option but it has some advantages also. Theoretical description 
is simpler for pure rotational excitation. The levels positions are better known for ground 
vibrational states. From the experimental side, it is easier to find a radiation having proper 
frequency because of better frequency tunability of microwave sources. 

We start with brief qualitative description of the phenomenon. CH 3 F has two types of 
states, ortho and para, shown in Fig. [I]. The data in this figure correspond to the 13 CH 3 F 
molecule and have been calculated using the molecular parameters from [fT4| |. Spin isomers 



of CH 3 F are distinguished by the total spin of the three hydrogen nuclei, 1 = 3/2 for 
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ortho and 1 = 1/2 for para isomers. For ortho isomers only rotational quantum numbers 
K = 0,3,6.. are allowed (K is the projection of molecular angular momentum, J, on the 
molecular symmetry axis.) For para molecules only K = 1,2,4,5... are allowed [H. 

There are two close pairs of ortho and para states in the ground vibrational state of 
13 CH 3 F that are significantly mixed by the intramolecular perturbation, V, and that are 
important for the ortho-para conversion in the molecule. For a qualitative description, let 
us take into account only one of these pairs, m — n, and assume that there is no external 
radiation yet. Suppose that the test molecule is placed into the ortho subspace. Due to the 
rotational relaxation caused by collisions, the molecule starts to shuttle up and down along 
the ladder of rotational states. Nonmagnetic collisions do not change the nuclear spin state 
directly, i.e., the relevant cross-section is zero, a(ortho\para) = 0. This shuttling along the 
rotational states inside the ortho subspace continues until the molecule jumps to the state 
m. During the free flight after that collision the intramolecular perturbation, V, admixes 
the para states n to the ortho state m. Consequently, the next collision has a probability 
(usually very small) to transfer the molecule to other para states. This localizes the molecule 
inside the para subspace and the spin conversion occurs. This is the mechanism of radiation 
free nuclear spin conversion induced by the intramolecular state mixing Jl5| (see also |16|). 

In case of a strong microwave radiation applied to the molecular transition q — n in the 
para subspace, mixing of the states is affected by the radiation which allows to control 
the ortho-para conversion. Influence of a radiation comes through the three major effects, 
level shift (dynamical Stark effect, well-known in nonlinear laser spectroscopy |L7| |I~8|j), level 
population change, and light-induced coherence. In general, these three components cannot 
be separated and work together. 

In order to consider a real molecule, the above simplified picture has to be developed 
further. One has to account the molecular center-of-mass motion. Although the intramolec- 
ular mixing does not depend on molecular velocity, the radiation-molecular interaction does. 
Consequently, the ortho-para state mixing in coherent control depends on molecular velocity 
too. 

Another complication comes from the degeneracy of molecular states. Even for the sim- 
plest case of pure radiation polarization (linear or circular) there are many excitation chan- 
nels each having its own absorption coefficient and saturation parameter. These channels 
differ by M-quantum number, projection of J on the laboratory axis of quantization. It is 
important also to keep in mind that there are other degeneracies of states. Each state of 
13 CH3F in Fig. [3] is determined by the set of rotational quantum numbers (J,K,M), total 
spin of three hydrogen, J, its projections on the laboratory z-axis, a, and ^-projects of spins 
of carbon and fluorine nuclei, both having spin 1/2. The energy of rotational states of CH3F 
depends only on J and K quantum numbers, if tiny hyperfine contribution to the level 
energy is neglected. We end the qualitative picture by summarizing important parameters 
of 13 CH 3 F in Table 1. 

Quantitative analysis of the isomer coherent control will be performed using kinetic equa- 
tion for the density matrix, p. The molecular Hamiltonian reads, 

H = -(h 2 /2m )V 2 r + H + hG + HV. (1) 

Here the first term is the Hamiltonian of the molecular center-of-mass motion with tuq being 
the molecular mass. The main part of the molecular internal Hamiltonian, H Q , has the eigen 
ortho and para states shown in Fig. [I]. hG describes the molecular interactions with the 
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external radiation that will be taken in the form of monochromatic travelling wave, 



G = - (E d/h) cos(uj L t - kr) , (2) 

where Eo, ujl and k are the amplitude, frequency and wave vector of the electromagnetic 
radiation, respectively; d is the operator of the molecular electric dipole moment. V is the 
intramolecular perturbation that mixes the ortho and para states in 13 CHsF. The mixing of 
m — n pair (7=9,^=3-11,1) is performed by the spin-spin interaction between the molecular 
nuclei |]TB|, |TP[, The pair m' — n' (20,3-21,1) is mixed by the spin-spin and spin-rotation inter- 
actions ]TB|, 120, 121], E2J. Account of the level degeneracy for the light-molecular interaction 



is a difficult problem, in general. It is most simple for the case of pure polarization, linear 
or circular. We will consider the electromagnetic radiation having linear polarization. 

In the representation of the eigen states of Hq (a-states) and classical description of the 
molecular center-of-mass motion, kinetic equation reads ||17||, 



dp/dt + v- Vp = S-i[G + V,p]. (3) 

Here S is the collision integral; v is the molecular center-of-mass velocity. Spontaneous 
decay is not included in this equation because it is negligible for rotational transitions in 
comparison with the collisional relaxation. 

Kinetic equation for the total concentration of para molecules can be obtained directly 
from Eq. (||), 

dp p /dt = -2Re J p m 'n'V n > m > + ^ pmnV nm )d\. (4) 

Here the total concentration of para molecules, p p = J2 a I Pp( a > v)<iv, a G para. Summation 
is made over all degenerate sublevels of the states m', n! and m, n. In Eq. (|]) a uniform spatial 
distribution of molecular density was assumed. Collision integral did not enter into Eq. (|J) 
because by assumption collisions do not change the molecular spin state, i.e., ^ a J S aa d~v — 
0, if a G ortho, or a G para. G did not enter into Eq. @ either because the matrix elements 
of G off-diagonal in nuclear spin states vanish. 

The off-diagonal matrix elements p mn and p m 'n' will be found in perturbation theory. 
Further we assume the perturbations V being small and consider zero- and first-order terms 
of the density matrix, 

P = P' + P". (5) 

Collisions in our system will be described by the model standard in the theory of light- 
molecule interaction. The off-diagonal elements of S have only decay terms, 

S aa ' = -TPaa'] Oi ^ a'. (6) 

The decoherence rates, T, are taken equal for all off-diagonal elements of collision inte- 
gral. This assumption simplifies the theoretical model. Note, that the dependence of the 
relaxation rates on rotational quantum numbers is known |23 . 



The diagonal elements of S are expressed through the kernel of collision integral, A, in a 
usual way, 

S(a, v) = J M a i v l a i> v i)p(«i> vi)dvi - p(a, v) ^ / A(a u v x \a, v)dv 1 . (7) 
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We consider the model of strong collisions with the following collision kernel for para 
molecules 



A(a, v\a\, vi) = v r w p (a)d(v — vi) + ^5 aai /(v); a, ol\ € para, (8) 

and similar equation for ortho molecules. w p (a) in Eq. (|j) is the Boltzmann distribution of 
rotational state populations of para molecules, 

w p (a) = Z' 1 exp(-E a /k B T), (9) 

with Z p being the rotational partition function; E a the rotational energy of a-state; T the 
gas temperature; ks the Boltzmann constant. The symmetry of CH3F is such that the 
partition functions for ortho and para molecules are practically equal at room temperature. 
Partition functions account all degeneracies of states (see Ref. [0 for more details); /(v) 
in Eq. (P) is the Maxwell distribution, 



/(v) = tt-^V exp(-v7^); ^0 = pk B T/m . (10) 

In Eq. (H), two relaxation rates were introduced, rotational relaxation (u r ) that does not 
affect molecular velocity and translational relaxation (u t ) that equilibrates velocity but do 
not change rotational state. Note that the rotational relaxation is accompanied in our model 
by the relaxation in M quantum numbers. Note also, that the collisions in the model do not 
change the molecular spin state. The introduction of different relaxation rates for different 
degrees of freedom makes the model of strong collisions more accurate and flexible. It allows 
to adjust the model to the particular experimental conditions. Because of its simplicity 
model of strong collisions is often used in laser physics and nonlinear spectroscopy, see, e.g., 
fT7| , p4| , 25]. Numerical values for the collisional parameters T,u r , and v t will be determined 
later. 



III. MICROWAVE ABSORPTION 

For the zero-order term of the density matrix one has the following kinetic equation, 

dp'/dt + v • Vp' = S' - i[G, p'}. (11) 

Electromagnetic field interacts with para molecules only. Consequently, ortho molecules 
remain at equilibrium in the zero-order perturbation theory. For the level populations of 
ortho molecules one has, 

p' (a,v) = (N-p' p ) Wo (a)f(v), (12) 

where N is the total concentration of molecules. 

Eqs. and (JTTJ) allow to deduce an equation for the stationary populations of 

para molecules, 

(y r + v t )f/ p {a, v) = v r w p {a)p' p {w) + v t f(w)p' p (a) + p' p p[5 aq - S an ], (13) 
where the excitation probability, p, is defined as, 



r 2r I G q n I 

PpP ~ 



r 2 + (n - kv) 2 



Pp(n, v) - p' p (q, v)l . (14) 
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In Eq. fljjP the notations were introduced, 

p' P ( v )= Pp(oi) = p' p (at,v)dv. (15) 

a£para 

Eq. QUip is written in the rotating wave approximation. Nonzero matrix elements of G 
(electric field has linear polarization along z-axis) are given by, 

Gqn = G(M)e* (kr_ ™ ) ; G(M) = E 10 (d~io) qn /2h, (16) 

where Q = u>l — oj qn is the radiation frequency detuning from the absorption line center, 
u qn ; the bar over a symbol indicates a time-independent factor; E%q and d\o are spherical 
components of the electric field and electric dipole moment vectors, respectively The 
matrix elements of d 1Q reads U, 

|(^) 9n | 2 =|rf(M)| 2 = (2J 9 +l)(2J n + l)(_^ I J A d\ (17) 

where (:::) stands for the 3j-symbol and d is the permanent electric dipole moment of CH 3 F, 
d=1.86D[g|]. 

Solution of Eq. fl!3D has no difficulty and can be presented in the form, 

Pp(«> v ) = p' P w P ( a )f( v ) +p' P [( T 2 -n)pi/(v) + Txp\[8 aq - 5 an ], (18) 

where the relaxation times are t\ = (z/ r + z/ t ) _1 , r 2 = z/" 1 , andpi = / pdv. We have separated 
here the field-induced contributions nonequilibrium in a and v and nonequilibrium only in 
a. Solution ([18]) shows that radiation affects the population of only two states, q and n. 
This is the consequence of the accepted simple model of collisions. 
The excitation probability can be found from Eqs. ( p| ) and (|I8"D, 

0.5Aw 

Pi 



ti(kR) 1 + t 2 -ti' 

rV(v) g, 

p rj + (!i-kv) 2 i?' 1 ' 

where the difference of the Boltzmann factors is Aw = w p (n) — w p (q); the saturation pa- 
rameter, k, and saturation intensity, S sat , are 

S = s f^L, ; (2°) 



Ssat' sat tom\d{M)\* 

the homogeneous line width is Y B = Ty/1 + k, and 

^-7 r| + (fi-kv)2- (21) 

This integral can be expressed through the probability integral, but for numerical calcula- 
tions performed in this paper it is easier to calculate it straightforward. 
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In a similar way one can obtain from the kinetic equation (|ll|) the off-diagonal density 
matrix element, 



^ = _< ^# [r + i(n ~ kv)1 - (22) 

We can adjust now parameters of the collision kernel @. Kinetic equation ([11]) describes 
a diffusion process with the diffusion coefficient, D = v\j1vt. Diffusion coefficient for CH3F 
is equal to D ~ 10 2 cm 2 /s at the pressure 1 Torr. This determines the velocity equilibration 
rate, v t = 4.4 • I0 7 s _1 /Torr. 

Attenuation of the radiation is given by hu L p p J2M Pi- Consequently the absorption 
coefficient, x(f2), is determined by the expression, 

x(n) = hu LPp s- 1 Y / Pi- (23) 

M 

In the low field limit (k — > 0) it is reduced to, 

Xlowity = ^-ppAwRiow^S^, (24) 

where Ri ow — h m s^o R is the Foigt profile of the absorption line. If T ^> kvo the absorption 
line is Lorentzian having the width equal T. Experimental data on x(^) f° r the rotational 
transition 11, 1 — > 12, 1 can be used to determine the value of T. Equally, broadening 
of any other rotational transition can be used because in our collision model (§) one has 
the same T for all off-diagonal density matrix elements. There are experimental results on 
the broadening of the ortho-para transition (9,3)— (11,1) obtained from the level-crossing 
resonances in 13 CH 3 F nuclear spin conversion. This experiment gave the value Y/P = 
1.9 • 10 8 s _1 /Torr |27, |28j that will be used in the present calculations. The last unknown 



P P V 



parameter, rotational relaxation, u r , can be determined, e.g., from the power saturation of 
the absorption coefficient. This information is not available and we assume v r = Y. This is 
reasonable, because the pressure broadening in molecules is determined mainly by the level 
population quenching, although this estimation for v r is probably too high. 

We can now demonstrate the model at work by considering the microwave absorption by 
13 CH 3 F. Absorption spectrum is determined by the selection rules J—*J+1,K—*K. The 
spectrum consists of groups of lines nearly equally separated by 50 GHz. Inside each group 
the lines, different in K, are rather dense. Spectrum near the line 11, 1 — * 12, 1 is shown in 
Fig. 0. The two spectra correspond to low radiation intensity and to S = 100 mW/cm 2 and 
the gas pressure equal 30 mTorr in both cases. 

Saturation intensity for the line 11,1 — > 12,1 is equal to 43 W/cm 2 (M = 11) and 
6.8 W/cm 2 (M = 0) at the gas pressure 1 Torr. S sa t is proportional to the pressure squared, 
thus at 30 mTorr, S sa t —6 mW/cm 2 . An example of the absorption coefficient saturation 
is given in Fig. |3|. Because the Doppler width of the transition is small, kvo = 0.74 MHz, 
low field absorption in the line center depends weakly on CH 3 F pressure if P ^ 100 mTorr. 
Another example of the saturation effect is shown in Fig. (|. Here the relative level population 
difference [p p [n) — p p (q)}/ p p Aw is given as a function of radiation intensity. One can see that 
radiation having S = 100 mW/cm 2 decreases the level population difference significantly. 
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IV. FIRST ORDER THEORY 



The kinetic equation for the first order term of the density matrix p" is obtained from 
Eq. (§, 

dp"/dt + v ■ Vp" = S" - i[G, p"\ - i[V, p'\. (25) 

Ortho-para conversion is determined by the terms p" m i n i and p" mn because zero order matrix 
elements off-diagonal in nuclear spins vanish (see Eq. (|j)). Radiation does not affect the 
levels m! and n'. Consequently, the matrix element p" m i n i is not different from the case of 
the field free conversion ||16|| , 

/Cn^^Wn'-vJ-PiKv)], (26) 
1 + too 

where u' = io m 'n'- The density matrix element p" mn can be obtained from the equations 
which are deduced from Eq. 



(<9/<% + v-V + Y)p" mn -ip" mq G qn = -iV mn [p' p {n, v) - p' {m, v)]; 
(d/dt + v • V + T)pl q - i P " mn G nq = -iV mn p' nq . (27) 

Substitutions, V mn = V mn e^, (u = uj mn ); p" mn = p'^; p" mq = p'^ n+ ^-^, 
transform Eqs. (E7I) to algebraic equations which can be easily solved. The density matrix 
element that one needs for the kinetic equation (Q) reads, 

_. V F + + U - kv )] Wjh V ) - Po( m , V )] + iGqnPnq /^x 

Pmn ~ mn (T + tuj)[T + i(n + uj-kv)] + \G(M)\i ■ [ ' 

Note that indices m, n, q in Eqs. ([28]) represent the set of quantum numbers that comprise 
all degenerate sublevels. Consequently, one has nonzero terms ~p" mn for the combination of 
quantum numbers which are allowed by the selection rules for V. 



V. CONVERSION RATES 



Solutions for p" m t n i and p" mn together with the level populations from Eq. (11 
off-diagonal matrix element from Eq. (|22|) allow to present Eq. (|j) as, 

dp p / dt = iV( 7 ' + 7po ) - p p7; j = j'+ V + 7po + 7op - 7„ - 7, 



and the 



(29) 



In this equation we have neglected the small difference between the total concentration of 
para molecules, p p , and its zero order approximation, p' p . The partial conversion rates in 
Eq. (|29|) have the following definition. The field free conversion rate through the upper level 



pair, m' — n', 



po 



E 



2riK r) /„/i 2 



-WoVm 



(30) 



r 2 + u 12 

Equation for 7^ p is obtained from 7 po by substitution the Boltzmann factor w p (n') instead 
of w (m'). Summation is made here over all degenerate substates of m' and n' states. The 
rate 7 po is given by, 



■ 2T 



to' 



+ ReJ FJ(v)dv 



w (m 



(31) 



S 



Equation for 7 op is obtained from 7 po by substitution the Boltzmann factor w p (n) instead of 
w (m). The rates 7 po and 7 op are field dependent. Their zero field limits coincide with the 
field free conversion rates through the pair of states m — n. The "non-coherent" contribution 
to the conversion, the rate j n , originated from the radiation-induced level population change 
in Eq. (EH), is given by, 



In 



2Tr 2Pl 
T 2 + uj 2 



+ 2Re J [(r 2 - n) Pl f(y) + r xP \F x dv 



(32) 



And finally, the "coherent" contribution to the conversion rate, j c , originated from the p; 
in Eq. (EH) is, 



nq 



Pi 



T 2 + u< 



Re J pF 2 dv 



= 12 

ic / j I v mn | 

In Eqs. (pT]),(|32]), and ( p3[) the following functions were introduced 



(33) 



Fx 
F 2 



r + iuj ) Y\ + i(Q + io\ — kv) ' 

r + T 1 + iux 1 

r + iuj T\ + i(Q + to\ — kv) 




\G(M)\ 2 

t 2 + uj 2 

\G(M)\ 2 \ 
T 2 + uj 2 J ' 



(34) 



Interpretation of the introduced conversion rates can be considered as follows. Strong 
resonant radiation splits the molecular states, change the levels populations and introduce 
coherences in the molecule. The field dependent part of 7 po can be considered as due to the 
radiation induced level crossing. This term has resonance at Q = —uj\. The first term of 7„ is 
due to the population effect. It has the resonance at Q = where the excitation probability 
has maximum. The second part of 7„ is due to the population change and level-crossing. 
It has two resonances, at Q = and at Q — —U\. The coherent contribution, 7 C , has also 
resonances at these two frequencies. 

The two peaks in the conversion rate spectra have rather distinctive features. The first 
one, at Q = 0, is quite similar to the radiation free conversion rate, e.g., 7 po (Eq. ([5t])). In 
our case u 3> T and the amplitudes of the peaks at Q = are proportional to T, thus to 
the gas pressure. In the limit to ^> I, contributions to the conversion rate provided by these 
terms are similar to the ordinary gas kinetic processes which are proportional to the gas 
pressure also. 

The resonances at Q = ui\ have completely different signature that would result from a 
field free conversion pattern of a degenerate ortho-para level pair. In this case the conversion 
rate has 1/T dependence, see, e.g., Eq. (|30|). It allows us to refer to the resonances at 
Q = uj\ as produced by the crossing the ortho and para levels and resulting from the applied 
electromagnetic field. To reveal this property of the resonances at ft — —UJ\, let us estimate 
the integral, 

/(v)dv 



Ft + z(Q + u x - kv) ' ^ 35 ^ 

at the radiation frequency = uo\. If the limit of large Doppler broadening is valid, T\ <C kvo, 
the integrand in Eq. (p5|) has sharp resonance at v = 0. One can substitute /(v) by /(0) in 
Eq. d35|) and obtain, I oc I/Uvq. Thus, the conversion is produced through the degenerate 
ortho-para level pair (level crossing) having the width equal the Doppler width, kvQ. In 
the opposite limit of large homogeneous broadening, T x ^> kvo, one can neglect kv in the 
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denominator of the integrand (35) and obtain / oc 1/ri. Again, this is the conversion 
through the crossed ortho and para states but having the width Ti now. 

It would be useful to compare the results of the present model with the qualitative model 
[ ITJ . This comparison cannot be made directly because in [|10[ rovibrational excitation of 
CH 3 F was considered. But one can apply the idea of light-induced enrichment solely through 
the level population change and make the comparison. In the present notations, the field 
effect from the level population change is given by the first term in Eq. 

12 

r 2 + ^2 



2T\V , 

' ran fnr\ 

r2 . , 2 TaPi- (36) 



It gives larger amplitude for the peak at Q = than the present model. In the present 
model one has partial cancellation of peaks at Q — 0, see the first terms in the expressions 
for 7„ and j c . Because of the assumption, T — v r , made in the collision model we have two 
times smaller peak at Q = in the present model. The cancellation is less significant if F is 
larger in comparison with v r . 

We turn now to numerical calculations of the conversion rates in 13 CH 3 F. Contribu- 
tion from w! — n' pair is difficult to calculate using Eq. ( |3~0|) directly because of some 
uncertainty in the parameters involved. Instead, one can use the experimental value, 
7p D = 2.3 ■ 1CT 3 s~ 1 /Torr |27, |28|] and scale it linear in pressure. Such pressure dependence 



for 7' is valid if T <C uo' (see Eq. (p0|) ) which is fulfilled for the pressures P < 10 Torr. 

Calculation of the rates 7 po , 7 n , and 7 C needs the matrix elements V mn . Mixing of the 
states m and n is performed by the intramolecular spin-spin interaction between the nuclei 
of 13 CH 3 F. Dependence of V mn on nuclear spin variables is accounted easily by summation 
because other factors in Eqs. (|3~T|)-(|3"3|) do not depend on nuclear spins. Then, the only 



remaining degeneracy is in M-quantum numbers. This quantity reads, 

]T \V mn \ 2 = (2J m + l)(2J n + 1) ^ q kI) { -M' M'-M M ) T ^ (37) 

Here 72 i9 (q = K m — K n = 2) is the magnitude of the spin-spin interaction and summation 
is made in all nuclear spin projections. The value of T2 2 calculated from the molecular 



structure is equal to 69.2 kHz. This value is confirmed by the experiment p9]| . But in 
this paper we have to take somewhat smaller value, 7^ 2 = 64.1 kHz, as it was obtained 
self-consistently for the three parameters, T 2i2 , T, and ^' po P8| . 

Examples of the conversion rates are shown in Fig. [5[ The upper panel gives the rate 7 
for two pressures, 30 mTorr and 100 mTorr. The radiation intensity in both cases is equal 
to 100 mW/cm 2 . One can see that the peak at Q = —uo\ is ~3 times larger at 30 mTorr 
than at 100 mTorr. There is also peak at Q = having "negative amplitude" but it is too 
small to be visible in the upper panel. The lower panel shows the field dependent rates, 
— {in + 7c)- They are taken with the same sign as they contribute to 7 in Eq. (^). One can 
see from this panel that the pressure dependences of the amplitudes of these two peaks are 
opposite. This is the consequence of the crossing of ortho and para states at fl = —uj\ and 
off-resonant nature of the peak at Q = 0. 

Broadening of the two peaks in the conversion rate are very different too. The peak at 
Q = has the broadening as ordinary absorption line does. At large saturation parameter, k, 
its width is ~ 2|G| and grows rather fast with intensity. The width of the peak at Q = —uji 
is given by ~ T(l + \G\ 2 /uj 2 ). Consequently, the power broadening of this peak is very small 
at our conditions. Fig. E illustrates the difference in the peaks widths. 
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VI. ENRICHMENT 



Solution of the kinetic equation (|29|) can be presented as, 

P P (t) =P P + (p P (0) - p p ) exp (-7*); p p = N(i po + 7 po )/7, (38) 

where p p and p p (0) are the steady-state and initial (equilibrium) concentrations of para 
molecules, respectively. Enrichment of para molecules will be defined as, 

««) = ^-i- (») 

Partition functions for ortho and para isomers of CH 3 F are equal. Consequently, 7 po = j' , 
l P o = 7o P and enrichment, /3(fi), can be expressed as, 

/3(fi) = ( 7 „ + 7c)/7, (40) 

which will be used for the numerical calculations of /3(Q). One can note from Eq. ([!(]) that 
despite the rates 7 po and 7 op have rather large field dependent parts (see Fig. (|)), they alone 
would not produce an enrichment. It can be understood because these field dependent parts 
are due to the mixing of states shifted by radiation but having equilibrium populations (see 
Eq. flSiP). Conversion through the mixing of equilibrium populated states does not affect the 
ortho-to-para ratio. Examples of the enrichment /3(Q) are given in Fig. || At the pressure 
30 mTorr and radiation intensity S = 100 mW/cm 2 , enrichment peak at Q = —uj\ is ~3% 
and it is ~4 times higher than the enrichment at Q = 0. At larger radiation intensity, 
amplitude of the enrichment peak, j3(—ui), saturates at ~5%. The peak at Q = grows to 
~1%. 

The simplified model that accounts only the radiation-induced level population change 
predicts a factor two higher peak at Q = than the present model. As was discussed above, 
in the present model one has partial cancellation of contributions originated from the rates 
7„ and 7 C which results in a smaller enrichment. Another significant difference between these 
models is that the simplified model does not predict an enrichment at Q — 

Broadening of the two peaks in enrichment is very distinctive and is similar to the peaks 
of conversion rate, although there is some difference. The enrichment peak at Q = —uji 
is broader than the peak of 7 because of large 7 po and 7 op at resonant frequency in the 
denominator of Eq. ( flO|) . 



VII. DISCUSSION AND CONCLUSIONS 

We have developed a model of spin isomer coherent control that accounts for the molec- 
ular level degeneracy in magnetic quantum numbers. This degeneracy together with the 
account of molecular center-of-mass motion allow now to apply the model to real molecules, 
having actual level structure and ortho-para mixing Hamiltonian. The developed model was 
used to analyse the microwave induced enrichment of spin isomers in 13 CH3F. Spectrum of 
the conversion rate consists of two peaks with the peak at Q = —uo\ being two orders of 
magnitude higher than the peak at Q = 0. 

The enrichment spectrum has also two peaks at the same frequencies. One can obtain 3% 
enrichment using forbidden (for ordinary absorption) resonance at Q = —U)\, gas pressure 
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30 mTorr at room temperature and radiation intensity S = 100 mW/cm 2 . Amplitude of the 
peak at Q = is 4 times smaller. Cooling the gas to 200 K would increase enrichment to 
~5% because of the increase of the level population difference. At higher radiation intensity 
enrichment saturates at ~5% if the gas temperature T = 295 K and at ~7% if T = 200 K. 

For the practical implementation of the microwave enrichment, the use of the peak at 
Q = —uj\ has a few advantages. Enrichment here is significantly larger than at Q = 0. 
Moreover, spurious effects can decrease the enrichment at Q = even further. It can be due 
to the gas heating by radiation. Note, that for the peak at Q — —oj\ absorption is negligible 
and there is no gas heating. Another spurious affect may be due to the resonant exchange 
of rotational quanta between ortho and para isomers. This effect prevents depopulation of 
one rotational state by radiation (state n in our case) in comparison with the population of 
the state m having nearly equal energy. Again, the peak at Q = —u)\ has the advantage of 
very low absorption coefficient and thus low radiation induced population change. Finally, 
the disadvantage of the peak at fl = is the spurious absorption by the line 11, — > 12, 
(Fig. H). Although at low pressures this line is well separated from the line 11, 1 — > 12, 1 
(the gap is 10 MHz), power broadening partially overlap these lines. On other hand, the 
peak Q = —uo\ is situated at the blue side of the line 11,0 — > 12,0 being 120 MHz away 
from the nearest absorption line (Fig. 0). 

Enrichment obtained by microwave excitation of 13 CH3F at reasonable experimental con- 
ditions is not large, ~3%. On the other hand, there are some applications where such 
enrichment could be significant, e.g., in spin isomer enhanced NMR technique |7|, M. Note 
that for the standard 200 MHz NMR the difference in Boltzmann factors between Zeeman 
states, which determines the amplitude of NMR signal, is only 3 • 10" 5 . If even a fraction of 
the 3% isomer enrichment would be transported to the Zeeman level populations it would 
enhance the NMR signal significantly. 

Apart from any possible applications of microwave enrichment, which is by far too early 
to discuss now, observation of the microwave enrichment would have prove-of-principle im- 
portance for the coherent control of spin isomers. It is interesting to note also that the 
isomer enrichment at Q = —uji would demonstrate an example of enhanced access to weak 



processes in molecules through the isomer enrichment ||30|| . Suppose, one would like to 
measure the absorption at Q — —u>i directly. At the conditions considered in the paper, 
x(—ooi) ~ 6 ■ 10~ 6 cm -1 which is rather difficult to measure. This value has to be compared 
with the 3% enrichment in the coherent control which should not be difficult to measure. 
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Table 1. Positions of levels in 13 CH3F. Molecular parameters are from Ref. 



Notation* 


J,K 


I 


E(cm- 1 ) 


Frequency (MHz) 


m' 


20,3 


3/2 


387.1 


351.01±0.16 {m! -n') 


ri 


21,1 


1/2 


387.1 




Q 


12,1 


1/2 


133.7 


596294.285±0.013 (q - n) 


n 


11,1 


1/2 


113.8 


130.99±0.15 (n-m) 


m 


9,3 


3/2 


113.8 


596425.28±0.15 (q - m) 



* Notation in Fig. [l|. 
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ortho 



para 



FIG. 1: Position of rotational levels of 13 CH3F. Molecular parameters are from Ref. 0]. Two 
pairs of states important for the spin conversion in this molecule are indicated. Small vertical line 
in the para subspace indicates microwave excitation of the transition n — > q. Two bent vertical 
lines indicate rotational relaxation. Parameters of the important states are summarized in Table 1. 
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Frequency (MHz) 



FIG. 2: Absorption spectrum near the line 11,1 — > 12, 1. Numbers in the graph indicate the 
K- values. The gas pressure is 30 mTorr. Low intensity absorption is shown by thin line, the case 
of S = 100 mW/cm 2 is shown by thick line. 
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FIG. 3: Saturation of the absorption coefficient in the line center of transition 11, 1 — > 12, 1. Gas 
temperature is T = 295 K. 
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Radiation intensity, S (mW/cm ) 



FIG. 4: Saturation of the level population difference normalized to the field free population differ- 
ence. Transition 11, 1 — ► 12, 1. Gas pressure, P = 30 mTorr; gas temperature, T = 295 K. 
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FIG. 5: Conversion rates. The upper panel shows the total conversion rate, 7, for the pressures 
P=30 mTorr (thin line) and P=100 mTorr (thick line). The lower panel shows the field dependent 
contribution, — (j n + 7 C ), at the pressures P=30 mTorr (thin line) and P=100 mTorr (thick line). 
The radiation intensity is S = 100 mW/cm 2 
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FIG. 6: Enrichment of para molecules, f3, as a function of radiation frequency detuning, Q,. Gas 
pressure, P = 30 mTorr (thin line); P = 100 mTorr (thick line). In both cases the radiation 
intensity S = 100 mW/cm 2 and gas temperature T = 295 K. 
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